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We show that the dependence of high-order harmonic generation (HHG) on the molecular orien- 
tation can be understood within a theoretical treatment that does not involve the strong field of the 
laser. The results for H2 show excellent agreement with time-dependent strong field calculations for 
model molecules, and this motivates a prediction for the orientation dependence of HHG from the 
N2 3<r 9 valence orbital. For both molecules, we find that the polarization of recombination photons 
is influenced by the molecular orientation. The variations are particularly pronounced for the N2 
valence orbital, which can be explained by the presence of atomic p-orbitals. 
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I. INTRODUCTION 

In this paper we establish a connection between 
photoionization/recombination (weak field processes) 
and high-order harmonic generation (a non-linear phe- 
nomenon in strong- field physics). The latter process is 
one of the most studied aspects of intense-laser physics 
because it serves as a source of coherent radiation at high 
frequencies 0, Q- High-order harmonic generation can 
be explained by a recollision mechanism Q: Close to 
the maximum of the electric field of a femtosecond op- 
tical laser pulse a molecule is ionized. A free electron 
wave packet enters the continuum and follows the elec- 
tric field of the laser. If the laser is linearly polarized 
the electron will approach the molecule again. The most 
energetic recollisions take place near the second zero of 
the laser electric field after electron release Q. Hence, 
the laser field at the time of recollision can be considered 
as small. The optical laser drives the electronic wave 
packet far away from the molecule (as compared to the 
size of the molecule) . Due to rapid wave packet spreading 
the electronic wave packet will recollide approximately 
as a plane wave with its momentum parallel to the laser 
polarization. Possible consequences of recollision are re- 
combination, elastic scattering or double ionization. In 
the recombination process a photon is emitted, prefer- 
ably parallel to the beam axis of the incident laser and 
with a frequency that is a multiple of the incident laser 
frequency, therefore called high harmonic. By appropri- 
ate superposition of different harmonics one can create 
attosecond pulses which may be utilized to probe fast 
atomic and molecular processes [J. If the return time of 
the electronic wave packet is well defined one can even 
think of using HHG itself as a probe for time-dependent 
processes. The ionization by the laser would represent 
the pump pulse and the recolliding wave packet would 
represent the probe pulse. The time between those two 
events is shorter than an optical cycle of the laser. This 
may open the door to the time-resolved investigation of 
very fast atomic and molecular processes, cf. also the 
method described in Ref. 0. 

In recent years there has been growing interest in HHG 
from molecules. The dependence on molecular orienta- 



tion has been studied experimentally Q and theoretically 
0. Considering the complexity of this process, theoret- 
ical investigations have been carried out mostly for H2 
and H^" until now. How to overcome this? As indicated 
above, at the time of recollision, when the radiative re- 
combination occurs, the electric field of the laser can be 
considered to be small for the highest harmonics. In the 
following we will use an approximation in which the influ- 
ence of the laser field on the recombination is considered 
to be even negligible so that the computational methods 
developed in the context of photoionization can be used. 
Although this will not cover all the dynamics of HHG, it 
should explain quite well dependencies of the high har- 
monics on the molecular geometry and orientation. As 
we will show, this is indeed the case. As a consequence 
it will be possible to describe HHG in much more com- 
plicated targets in the future by shifting the focus from 
the exact treatment of the time evolution towards the 
exact treatment of the final molecular interaction, the 
recombination in high-order harmonic generation. 



II. METHOD 

In the recombination process the electron approaches 
the molecular core and a photon is ejected, leaving the 
molecule predominantly in its ground state. The dy- 
namics of this process is contained in its transition am- 
plitude. Since recombination is (microscopically) time 
reversed photoionization one can use as recombination 
transition amplitudes the complex conjugated photoion- 
ization transition amplitudes. Furthermore, since we are 
here not interested in near threshold behavior we can cal- 
culate those photoionization transition amplitudes eas- 
ily using the frozen core Hartree Fock (FCHF) method : 
The molecular ground state wave function is derived in a 
self-consistent-field approximation. The state of the ion- 
ized molecule is then obtained by removing one electronic 
charge out of the orbital that is ionized. The molecule 
is not allowed to relax ('frozen core'). The photoelectron 
orbitals were obtained using an iterative procedure to 
solve the Lippmann-Schwinger equation associated with 
the one-electron Schrodinger equation that these orbitals 
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satisfy (for further details see 0). 

For the mathematical description of the process the 
density matrix formalism [3] has been applied (see also 
To be able do so we must model the recombina- 
tion process. We know the ground states of the neutral 
molecule and of the singly charged molecular core, |Ao) 
and |A C ), respectively. We assume that the molecular 
orientation m does not change during the process. We 



also know the state 



,(+: 



of the incident electron. (In 



the following we will sum over not resolved molecular vi- 
brational states and not resolved spin polarization states 
of the electron.) Naturally, the photon, with a frequency 
Lu p , will be polarized. Therefore the full density matrix p 
of the state after recombination in dipole approximation 
reads 

p = |m A eujp) (mA eu p \ 

= T mA c kW\/mA c kW T+, (1) 

where e is the polarization vector of the photon and T 
the transition operator, i.e., the dipole operator. The 
photon properties will be measured in a detector in a 
direction n p . For a perfect detector one gets by pro- 
jecting on the different polarization states, which are 
|A) e{| — 1) , |0) , |1)} in an arbitrary reference frame, for 
the matrix elements p(A, A') of the density matrix p 



p(X, A') = (mAonpAwp|p|mAon p A'wp) 



(2) 



A common description of photon polarization employs 
the Stokes parameters. The Stokes parameters are de- 
fined in a reference frame which z-axis is parallel to the 
photon momentum In this frame the right (left) cir- 
cularly polarized photon state is | + 1) (| — 1)). (|0) does 
not exist in this reference frame due to the transversal na- 
ture of the light.) The four Stokes parameters are : the 
total intensity I, the degree of circular polarization p$ = 
(i+i — I-x)/I and the two degrees of linear polarization 
Pi = (7(0°) - 1(90°))// and p 2 = (7(45°) - 7(135°))//. 
(<ys in I(ip) starts at the x-axis in the xy-plane.) In the 
reference frame of the Stokes parameters one gets 



I = p(l,l) + p(-l,-l) (3a) 

P3 = (p(i,i)- P (-i,-i)). r 1 (3b) 

pi = -{ P {i,-i) + p{-iA))-r x (3c) 

p 2 = -i Goa,-!)-^-!,!))-/- 1 . (3d) 



The electronic wave function ki^y can be expanded 

into spherical harmonics However, due to the non- 
spherical molecular potential the dipole selection rules 
do not restrict the expansion as in atoms. Nevertheless, 
both, bound and continuum electron wave functions con- 
verge quite rapidly. Therefore, to a very good approxi- 
mation a limited number of terms is sufficient, truncating 
the expansion at a certain l max . p(X, A') can be split into 
a kinematical and a geometrical part 
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where (., .|.) are Clebsch-Gordan coefficients. In Equ. 0]we have used l m ax = 10. The geometrical dependencies are 
expressed in bipolar spherical harmonics, 

d £ 

y r d f 7 (m, k°) = £ £ (dS, CM\T - 7 ) Y dS (m) Y CM (k°) , (5) 

8=-dM = -C 

with spherical harmonics Yk q and k° as the normalized electron momentum. The reference frame is given through 
the photon. The dynamical coefficient is 

M dCr = dtY, A £ EE I) E H 1+mm+A ™ +r II' Rlm m X m R*>m' m \> m 

1=0 m m =-(A° + l) A m = -1 l'=0m' m = -{A° + l) X' m =-l 

( 6 ) 
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where A is if the recombined orbital has a symme- 
try and 1 if it has it symmetry. In Eq. JSJ Wigner 3J 
symbols have been used. A hat over a quantum number 
means I = \/2l + l. The dynamical part is calculated 



in the molecular body frame (symbolized by a sub-index 
m at the quantum numbers). By applying microscopic 
time reversal, the recombination matrix element in the 
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molecular body frame in length form is 

Rim m \ m = w p exp(iA/) (A c lm m \ d\ m |A A m ))* , 

(7) 

where lu p is the energy of the photon, A; is the Coulomb 
phase shift, and (A c lm m \ d\ m |AoA m ) is the photoioniza- 
tion dipole matrix element in the body frame with the 
dipole operator d\ m . 
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III. RESULTS FOR H 2 

Using the formulae of the last section one can calculate 
the photon intensity and polarization, the Stokes param- 
eters, as a function of the electron energy and for differ- 
ent orientations of the molecule and electron and photon 
propagation directions. Here we will focus on the HHG 
geometry (see Fig. where the electron momentum is 
perpendicular to the emission direction of the photon. In 
the following, we distinguish between two cases: (I) the 
molecule lies in the plane spanned by the electron mo- 
mentum and the photon direction, and (II) the molecule 
rotates in the plane perpendicular to the photon direc- 
tion. We have first calculated the intensity for H2 recom- 
bination in geometry (I) as a function of the molecular 
orientation. Here, the photon is polarized parallel to the 
electron momentum for symmetry reasons. The calcula- 
tion was carried out for different bond lengths and for 
different wavelengths of the electron (see Fig. |2J). One 
finds a pronounced minimum, which shifts if one changes 
the bond length of the molecule or if one changes the 
wavelength of the electron. 
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FIG. 1: The recombination geometry is shown schematically. 
The molecule with a bond length ro is oriented relative to 
the electron momentum k e under an angle 9. The photon is 
emitted perpendicular to the electron momentum. 

One can explain the general behavior at electron wave- 
lengths comparable to the bond length of the molecule by 
the well-known two-center interference. Here one imag- 
ines the diatomic molecule as two centers which are hit 
coherently by the same plane electron wave, but with a 
phase difference that depends on the molecular orienta- 
tion towards the electron. Recombination leads to the 
ejection of a photon. Interference occurs, since it is not 



FIG. 2: Dependence of the recombination photon intensity 
on the angle between molecular orientation and the electron 
momentum for H2. The molecule lies in the plane spanned by 
electron momentum and photon direction [geometry (I)] . The 
solid curves are for an electron wavelength A e of 1.4 a.u., the 
dotted curves for \ e = 1.6 a.u. and the dashed curves for A e = 
1.8 a.u. (a) Molecular bond length ro = 1.4 a.u.; (b) ro = 
2.0 a.u.. The positions of minima and maxima are marked by 
arrows. Obviously, the positions of the extrema depend on 
the bond length and on the wavelength of the electron. This 
behavior can be explained by a two-center interference model 
(see text). Clearly, the minima are very pronounced. 

known which center has emitted the photon. Changing 
the electron wavelength and/or the molecular orientation 
will alter the phase difference so that an interference pat- 
tern will be obtained. At the energies used here, the pho- 
ton wavelength is much larger than the dimension of the 
molecule, so that one can neglect the phase shift resulting 
from the orientation of the molecule with respect to the 
photon. The bond length ro of the molecule, the electron 
wavelength A e and the angle 9 ex between molecular axis 
and propagation direction of the electron under which ex- 
trema in the recombination photon intensity appear (see 
Fig. ^) are then - in the two-center interference picture 
- related through 

TL 

r cos(8 ex ) = - A e + Sg, n = 0,l,2,... (8) 

where Sg is the difference of additional phase shifts the 
electronic wave function experiences in the vicinity of the 
nuclei. In the ideal case those phase shifts are equal and 
Sg is zero. 5g depends on the orientation of the molecule 
and is expected to be maximal if the molecule is parallel 
to the electron momentum and zero if perpendicular. If 
Sg is small, interference will be constructive for even n in 
Eq. © and destructive for odd n. Parallel and perpen- 
dicular orientation of the molecule relative to the electron 
momentum always give rise to trivial extrema. For fixed 
bond length ro and increasing electron wavelength, min- 
ima will occur at positions where the molecule is more 
and more aligned along the electron momentum, up to 
the point where both are parallel. In the process the 
minimum gets less pronounced and its absolute value is 
not approximately zero anymore. 

A convenient way of analyzing the extrema is presented 
in Fig. |2| where the projection r cos(9 ex ) is plotted as a 
function of A e . Our results bear a strong resemblance 
to those of time-dependent strong field calculations of 
HHG in H2 and H^ model molecules [l2|. This supports 
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FIG. 3: Recombination in geometry (I) for H2. For the 
extrema in the photon-intensity orientation dependence, the 
projection of the bond length on the electron momentum di- 
rection, ro cos(6ex), is plotted against the electron wavelength 
X e . The dotted (dashed) curve marks minima (maxima) due 
to two-centre emittance (see text). The solid curve represents 
minima for ro = 1.4 a.u., whereas the long-dashed (dashed- 
dotted) curves represent minima (maxima) for ro — 2 a.u.. 
For comparison the extrema in the orientation dependence in 
HHG for time-dependent strong field calculations for model 
molecules [l2T | are plotted as well: interference minima for ( V) 
at r = 2 a.u., laser intensity I = 10 15 W/cm 2 ; (•) H 4 " at 
r = 2 a.u., I = 5 X 10 14 W/cm 2 ; (□) H 2 at r = 1.4 a.u., 
I = 5 x 10 14 W/cm 2 ; (+) interference maxima for H^" at 
r = 2 a.u., I = 10 15 W/cm 2 . Note that the present defi- 
nition of A e differs from [la . 



our prior assumption on that one can treat the recom- 
bination in HHG approximately as a weak-field process. 
Also, one finds the predictions made about the increase 
of Sg towards parallel molecular orientation confirmed. 
Not surprisingly, the signature of two-center interference 
fades with increasing electron wavelength. However, at 
the electron wavelengths considered, this effect can be at- 
tributed mainly to the decreasing ratio of kinetic energy 
of the electron to the ionization threshold. In general, the 
orientation dependence of the recombination photon in- 
tensity for H2 can be well described within the two-center 
interference model. 

As indicated, we can calculate all the Stokes param- 
eters. In geometry (I), only linear photon polarization 
is possible. In geometry (II), however, where the molec- 
ular axis lies in the plane perpendicular to the photon 
propagation direction, the photon can have different po- 
larizations and even circular polarization can be obtained 
(see Fig. @J. All polarizations show strong variations in 
the vicinity of the interference minimum. Otherwise only 
small polarization variations have been found. Thus, al- 
though the polarization depends on the geometry, the 
difference is small for H2 because the signal is dominated 
by the polarization parallel to the electron momentum 
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FIG. 4: The H2 molecule lies in the plane perpendicular to the 
photon propagation direction which is parallel to the z-axis 
in the Stokes parameters frame [geometry (II)]. The electron 
is here chosen to move along the x-axis. Shown is as solid 
curve the recombination photon intensity for A e = 1.6 a.u., 
ro = 1.4 a.u. against the angle between the molecular axis and 
the propagation direction of the electron. Dotted, dashed- 
dotted and dashed curves are the Stokes parameters pi, P2 and 
pz, respectively. These show strong variations in the vicinity 
of the minima. 

except in the small range around the interference min- 
imum. This can be understood within the two-center 
model since the H2 molecular orbital is approximately 
the sum of two atomic ls-orbitals. These are spherically 
symmetric and therefore do not produce a signal polar- 
ized perpendicular to the electron momentum. 



IV. RESULTS FOR N 2 

Given the excellent agreement of our H2 results with 
time-dependent strong field calculations [L2| we can move 
to a prediction for the orientation dependence of HHG 
from the N2 3a g valence orbital. The time-dependent 
HHG calculation for this system is quite complicated and 
has not been carried out. 

While both, H2 la g and N2 3<r g , have the same sym- 
metry, they are rather different otherwise. While \a g is 
mainly built up from atomic s-orbitals and does not pos- 
sess nodes, 3<r g is dominated by atomic p-orbitals and 
has a more complex structure |l31 |. As a consequence 
the orientation dependence for N2 3a g is more complex 
than for H2 la s . As in the previous section, we have 
investigated geometries (I) and (II) . In Fig. [5] the ex- 
trema for the equilibrium bond length 2.068 a.u. as well 
as for 1.768 a.u. and 2.368 a.u. are plotted. Contrary to 
H2, there are big differences between the two geometries. 
Figure [HI shows the orientation dependence of the Stokes 
parameters for geometry (II) . Large variations are found 
over a broader range of angles than in H2, i.e., the com- 
ponent perpendicular to the electron momentum cannot 
be disregarded. At small angles, the signal is still dom- 
inated by the polarization parallel to the electron, but 
not so for larger angles. 
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FIG. 5: Projection ro cos(8 ex ) versus electron wavelength X e 
for molecular orientations 9 ex under which minima (x) and 
maxima (+) occur in the N2 3<t 9 recombination photon inten- 
sity. The bond length ro is 1.768 a.u. in (a) and (d), 2.068 a.u. 
in (b) and (e) (ground-state bond length), and 2.368 a.u. in 
(c) and (f), respectively. The upper plots [(a), (b), (c)] show 
the positions of extrema when the molecule lies in the plane 
spanned by the electron and photon direction [geometry (I)], 
whereas the molecule rotates in the plane perpendicular to the 
photon direction in the lower plots (d), (e), (f) [geometry (II)]. 
The dotted (dashed) curve is the one for minima (maxima) ac- 
cording to the two-center interference model for p-orbital con- 
tributions. Note that minima and maxima are interchanged as 
compared to the s-orbital contributions. The dashed-dotted 
line marks minima according to ro cos(8 EX ) = 0.75 A e . 



Since N2 is a homonuclear diatomic molecule we might 
expect to find two-center interference in the region where 
the wavelength of the electron equals approximately the 
internuclear distance. However, Fig. [S] shows that it is 
not straightforward to identify such signatures. 

To understand the observed behavior, we first note 
that the two atomic p-orbitals "inside" the N2 valence 
orbital are of course not spherically symmetric. There- 
fore, unlike s-orbitals, each of them can produce a sub- 
stantial component polarized perpendicular to the elec- 
tron momentum with a pronounced dependence on the 
molecular orientation. However, this component does not 
show up in geometry (I). This explains the difference 
between the two geometries. Furthermore, the molec- 
ular orbital is not constructed of atomic p-orbitals only, 
but we have an s-orbital admixture of about 30%. This 
makes the two-center interference picture problematic be- 
cause different interference patterns are expected for the 
two orbital types: to ensure a g symmetry, the two s- 
orbitals <fi s (r — ro/2) and (j) s (r + ro/2) are added with 
the same sign, so that interference conditions are ob- 
tained as described for H2; the p-orbitals 4> p (r — ro/2) 
and 4> p (r + ro/2), on the other hand, are added with op- 
posite signs, leading to an interchange of maxima and 
minima |12j . 

Clearly, the simultaneous presence of both types of in- 
terference will lead to a complicated orientation depen- 
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FIG. 6: The N2 molecule lies in the plane perpendicular to the 
photon propagation direction which is parallel to the z-axis in 
the Stokes parameters frame [geometry (II)]. The electron is 
here chosen to move along the x-axis. Shown is as solid curve 
the recombination photon intensity for A e = 1.6 a.u., ro = 
2.068 a.u. against the angle between the molecular axis and 
the propagation direction of the electron. Dotted, dashed- 
dotted and dashed curves are the Stokes parameters pi, P2 
and P3, respectively. 



dence. However, we can look for situations where either 
s- or p-orbitals dominate the signal. One such case is 
geometry (I) near an orientation of 90°. Here, the indi- 
vidual p-orbitals generate a negligible signal due to their 
mirror antisymmetry, as is obvious when one considers 
a matrix element of the form (0 Px (r)|z| exp(ifcz)) where 
the incoming electron is approximated by a plane wave. 
Consequently, s-orbital contributions should dominate. 
In fact, Figs. Ela)-(c) show that in the small-wavelength 
regime, geometry (I) yields local maxima at 9 = 90° 
as a consequence of constructive interference (similar to 
H2). Another example is geometry (II), when only the 
component perpendicular to the electron momentum is 
measured. In this case, the s-orbital contributions are 
small as explained above in the context of H2. Conse- 
quently, we should observe the p-type interference pat- 
tern, i.e., zero signal at 90° and a series of minima and 
maxima when the angle is decreased. This is indeed the 
case for small electron wavelengths as is shown in Fig. 
where the extrema for the perpendicular component 
are plotted. In this plot, the extrema are systematically 
slightly below the "perfect" two-center interference lines 
ro cos(8 ex ) — n A e /2. The total photon intensity in geom- 
etry (II) exhibits a similar behavior, see the lower panels 
of Fig. [S] This demonstrates the predominance of the 
p-orbital part even in the total signal. 

In geometry (I) we have both s- and p-orbital contri- 
butions for the small and intermediate angles. Although 
in this regime the results cannot be explained in a simpli- 
fied picture, we find a set of minima (see Fig. |SJ) following 
a straight line rocos(6 ex ) = 0.75 A e that lies just in the 
middle between the two-center interference lines. 
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FIG. 7: Same as Fig. 03 when only the component polar- 
ized perpendicular to the electron momentum is measured in 
geometry (II). (The molecular axis lies in a plane perpendic- 
ular to the the photon propagation direction.) The dotted 
(dashed) curve is the one for minima (maxima) according to 
the two-center interference model for p-orbital contributions. 

V. CONCLUSIONS 

In conclusion, we have shown that the orientation de- 
pendence of the recombination photon intensity in H2 



can be described very well in a two-center interference 
model. Our results on the orientation dependence bear 
a remarkable resemblance with those obtained from the 
time-dependent Schrodinger equation for HHG 01 • This 
shows that our method can be used alongside those oth- 
ers to obtain estimates about the effects of molecular ge- 
ometry and orientation on the photon intensity in HHG. 
We have made such a prediction for the case of N2 . Fur- 
thermore, we have demonstrated, that the photons from 
HHG in oriented molecules do not exhibit only linear 
polarization. Rather, the polarization of those photons 
show strong variations and can even be circular, depend- 
ing on the molecular orientation. For N2, the interpreta- 
tion of the results within a two-center interference picture 
is hampered by the fact that the valence orbital has ad- 
mixtures of both atomic s- and p-orbitals, which produce 
different interference patterns. However, we have pointed 
out situations where one of the two orbital types domi- 
nates the signal so that interference can be observed. 
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